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We analyze the coherent dynamics of excitons in three dimensional topologically disordered networks with 
traps. If the interactions between the nodes of the network are long ranged, i.e., algebraically decaying as a 
function of the distance between the nodes, the average survival probability of an exciton surprisingly shows a 
characteristic decay with features similar to the decay found for regular one-dimensional systems. We further 
show how this decay can be related to the eigenstates of the same system without a trap. 
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I. INTRODUCTION 

Coherent dynamical processes are of interest in many fields 
of research. Naturally, quantum mechanical systems like 
atoms in optical lattices display coherent dynamics, see QJJ] 
and references therein. Optical systems, for instance, lattices 
of evanescently coupled waveguides allow to study the basic 
principles of coherent energy transfer from the experimental 
side, given that these systems have long decoherence times 
10]. Also some biological processes show instances of being 
coherent: an example is the highly efficient energy transfer 
in photosynthesis, which cannot be understood based only on 
(incoherent) diffusion yl |4j, |5|] . 

The transport takes place over many, quite diverse geome- 
tries, ranging from very ordered lattices to random networks 
of active sites. Here, an efficient approximate description of 
the dynamics is given by a tight-binding ansatz which also al- 
lows to take the structure over which the excitation dynam- 
ics occurs fully into account la, tZl |8[]. Depending on the 
system, the interactions between sites extend from nearest- 
neighbor couplings to being of long-range type. An example 
for systems which allow to experimentally control and manip- 
ulate the range of the interactions are ultra-cold Rydberg gases 
lUlQli. Furthermore, it is possible to excite individual atoms 
locally; at later times these excitations may spread over the 
whole aggregate. By considering only the original Rydberg 
state and the state the local excitation is in, one can view the 
whole gas of atoms as being a collection of two-level systems 
ifTTll ; this resembles the situation in so-called spin-gases lfl2ll . 
The typical interactions between the atoms of the Rydberg gas 
is of dipole-dipole nature, decaying with the third power of the 
mutual distance between the atoms. 

Now, instead of directly following the propagation of a sin- 
gle excitation one can monitor the (global) transport process. 
One can introduce dissipation into the system in a controlled 
way. So-called excitation traps, which are placed at certain 
sites of the network, will "absorb'V'trap" the excitation once 
it hits the trap. For Rydberg atoms such traps can be realized 
by exciting one atom into a different state such that the exci- 
tation gets to be absorbed at this site lfl3l [Till . By monitoring 
the decay of the survival probability for the excitation to re- 
main inside the system, i.e., not to get absorbed by a trap, one 
obtains information about the dynamics within the system. 

The paper is organized as follows: After briefly recall- 



ing the phenomenological approach to excitation trapping, we 
show how the exciton survival probability for a given system 
can be related to the eigenstates of the same system without 
traps. Our model system, a random assembly of two-level 
systems in a cubic box, is described next, before presenting 
our numerical results, which were obtained using LAPACK 
routines. 



II. COHERENT EXCITON TRAPPING 

The Hamiltonian for a system of N nodes without traps is 
Ho, where the states | j) are associated with excitations local- 
ized at the nodes j (J = 1, . . . , N). These states are assumed 
to form a complete, orthonormal basis set of the whole acces- 
sible Hilbert space ((k\j) = Skj and ^ fe \k )( k\ — 1). Within 
a phenomenological approach, the Hamiltonian, which incor- 
porates trapping of excitations at the nodes m (m S M and 
M C {1, . . . , N}), is given by H = H iT, where iT = 
*r EmeM I to) (to I is the trapping operator, see Ref. lfl3ll for 
details. As a result, H is non-hermitian and has N complex 
eigenvalues, Ei — ei — vyi (I — 1, . . . , N) where 7; > 0, 
and N left and N right eigenstates, denoted by \*f?i) and 
respectively. The transition probabilities are given by 

Kkj(t) = |^exp(-7^)exp(-ie^)(fc|^)(# i |j)| 2 , (1) 
l 

so that the imaginary parts 7; of E\ determine the temporal 
decay. 

The mean survival probability Hm (t) of an excitation in the 
presence of M trap nodes is a global property of the network 
and is defined as 

= N^M £ £ (2) 

i.e., Hm (t) is the average of 7Tjy (t) over all initial nodes j and 
all final nodes k, neither of them being a trap node. For inter- 
mediate and long times and for a small number of trap nodes, 
njvf(i) is mainly a sum of exponentially decaying terms [13]: 

1 N 

Hm W « j^-jj cxp(-2 7 ;i)- (3) 
1=1 
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If the imaginary parts 7; obey a power-law with an exponent 
mean survival probability scales at in- 
termediate times as Hj^(t) ~ 

In the subsequent calculations we will use the following 
(approximate) average survival probability for a single trap 
(M = 1), obtained by slightly modifying Eq. {3), 

1 N 

n w^E cx p(- 2 ^)' w 

i=l 

which ensures the same normalization for different N, i.e., 
n(0) = 1 for all N. 

A. Relation to the eigenstates of Ho 



we consider interactions decaying as A~^. In the absense 
of traps the corresponding Hamiltonian H has the following 
matrix elements 

f-Ajl for/,- , j 

{k\U \j) = Ij^Ajl for k = j. (8) 
[kjtj 

We now choose one of the N nodes to be a trap, i.e., for this 
node the full Hamiltonian H has an additional purely imagi- 
nary matrix element —iT. Since the configuration of nodes is 
random, we can (without any loss of generality) assume in the 
following that the node labeled 1 is the trap. Figure Q] shows 
one realization of the system in which the trap node is explic- 
itly labeled. 



In a perturbation theoretical treatment, assuming the trap- 
ping strength T to be small, the imaginary parts 7; can be 
related to the eigenstates ) of the system without traps 
lfl5ll . For a single trap at node k one has 



11 =r](fc|* 



(0)\ 



(5) 



In the continuum limit the index I becomes a continuous 
variable, say, x and the imaginary parts lie in the interval ji £ 
[0, 7max]- If we know the density p[j{x)] — dx{"f)/di of the 
7(2;), Eq. (0]i can be recast into the form 



n(t) 



dl p{l) exp(-2 7 t) = £{p(j)}, 



(6) 




where £{,0(7)} denotes the Laplace transform of p(j). On 
the other hand, if we know H(t), we can obtain the density 
p{l) by the inverse Laplace transform C^ 1 {Ii{t)} . 

Take as an example a linear system with a trap at one end. 
At intermediate times, the mean survival probability will de- 
cay as U(t) ~ t- 1 / 2 1 15]. From Eq. © one immediately 
obtains p(~f) ~ 7" 1 / 2 and from this by integration 7 ~ x 2 . 
We will elaborate on the relation between trapping and the 
eigenstates of a closed system elsewhere in more detail (F" 



III. QUANTUM NETWORK 

We start from a random configuration of (N — 1) identical 
nodes and one trap node. All N nodes are placed at random in 
a 3-dimensional box with Cartesian coordinates {x^}, with 
i = 1,2, 3. Then the distance between two nodes j and k is 
given by 



1/2 



(7) 



where the coordinates xj and x^' are homogeneously dis- 
tributed random numbers in the interval [0, N]. To relate our 
results to the energy transfer dynamics within Rydberg gases 



FIG. 1: (Color online) Exemplary random configurations of N — 
100 nodes in a cube, one of which is a trap (labeled red dot). 

For disordered systems, we calculate averages over B dif- 
ferent realizations following 



1 R 

( ' )r = "ft [ ' ] r' 



(9) 



where [ • ] denotes the realization r. In our calculations we 
also assume the trapping strength [r] r to be realization de- 
pendent, because we require it to be proportional to the diago- 
nal element [(1 |Hn |1)] in that particular realization, namely 
[r] = r[(l|H |l)] . As it will turn out, the dependence of 
the decay on the value of T is quite weak - different T mainly 
rescale the time axis. Thus we will consider only the two ex- 
treme cases: (a) T = 10~ 6 , for which a perturbation theoret- 
ical treatment can be justified, and (b) V = 1, such that the 
average trapping strength is of the same order as the diagonal 
elements of Ho at the node of the trap. 

While each realization [ll(t)] r corresponds to a specific 
spectrum of the [7;] , the correspondence of to the 

average (7;) 7? is not that straightforward. However, for all 
t, the function exp(— 27^) is convex, therefore, Jensen's in- 
equality applies, see paragraph 12.41 of 11711 . such that we 



3 



obtain for a given I 

< exp(-2 7l i)) fl > exp ( - 2i( 7i ) fl ). 
From this we get a lower bound for (ll(t)} R : 

N 

X 



(10) 



1 N 

(n(t)) iJ >-^exp(-2t( 7i )fl). (11) 



i=i 



IV. NUMERICAL RESULTS 



For a linear system with equal distance between neighbor- 
ing nodes it has been shown in Refs. 1131 [1511 . that the decay 
of Hm (t) shows characteristic features of quantum transport, 
such as an algebraic decay of the survival probability at in- 
termediate times. It is not a priori clear that these features 
will persist in higher dimensions or in disordered systems with 
long-range interactions. As we proceed to show, the quantum 
network for a disorded system in three dimension, as defined 
in the previous section, still displays, in the ensemble average, 
distinct temporal regions with different behaviors of the decay 
of(U(t)) R . 




FIG. 2: (Color online) Decay of [n(i)] as a function of rescaled 
time tT/N 3 for R = 500 different realization and N = 100: (a) 
for r = 10~ 6 and (b) for V — 1. The dashed red curves show the 
ensemble averages (II (£)}#. 

As an example, Fig.|2]shows [ll(i)] r obtained from Eq. (0]i 
for N = 100 and for i? = 500 different realizations. Fig.[2ja) 
gives the results for r = 10~ 6 andFig.|2|b)forr = 1. In both 
figures the average (U(t))n is shown as a dashed red line. Ob- 
viously, the various realizations differ strongly in their behav- 
ior, displaying stronger and weaker decays. However, in the 



ensemble average, the decay of (II(t))_R is qualitatively com- 
parable to the one found for regular (one-dimensional) sys- 
tems, in that it follows a power-law at intermediate times. Ad- 
ditionally, all realizations show monotonic decays, although 
with variations over several orders of magnitude in time. 
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FIG. 3: (Color online) Ensemble averages (II(t)) r for different N: 
(a) for T = 10~ 6 and (b) for V = 1. The scalings in the intermediate 
time regions for N = 100 and N = 1000 are shown as solid black 
lines along with the appropriate scaling law. The arrows are guides 
to the eye pointing at the bend of (n.(t)) R for N = 1000. 

An analysis of (IL(t))n is shown in Fig. [3] at the two val- 
ues r = 10~ 6 and r = 1 for different N. In all cases the 
intermediate time decay can be fitted by a power-law 

(n(*)) H ~ (12) 

where, different from the regular linear case lfl3l,[l5ll . the ex- 
ponent r) is now iV-dependent. It turns out that approximating 
the exponent by 



r,(N) = TjoiV", 



(13) 



while keeping T fixed, yields satisfactory results. Note, how- 
ever, that r/o and also /j, might still depend on T. 

Now, we can estimate fi based on the results of rj for N% = 
100 and N 2 = 1000. From tj = T)(Ni)/N? = r}(N 2 )/Ng 
one obtains 



m^A^-ln^jVx) 
In 7V 2 -hxNi 



(14) 



From the numerical values given in Fig. [3] we get approxi- 
mately n -0.166 and tj = 0.0349 for T = 10~ 6 and 
H Pa -0.204 and r; = 0.0313 for T = 1. 

From Fig. [3] as well as from Eq. ( fT3l we see that the ex- 
ponent T)(N) decreases with increasing N. Certainly, if N 
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becomes very large it is quite improbable (in the ensemble av- 
erage) for an exciton to encounter the single trap. Therefore, 
the decay of (H(t))u can only be observed at very long times. 

A. Relation to the eigenstates | °' ) 

The fact that there is a decay at all can be explained by 
recalling that for small T, the decay rates 7; are related to 
the eigenstates ) of the system without traps, see above, 
which in the ensemble average we can write 

< 7i } fl , = (r|(fc|*(°>}| 2 ^. (15) 

Now, imagine there are states l^j ^) which have no overlap 
with the trap node at position k. This would cause vanishing 
7; for some I, which would lead to a non-vanishing H(t) even 
in the asymptotic time limit, see Eq. @. From Fig. [2] one 
sees that at least in the displayed time interval all r de- 

cay strictly monotonically. Therefore, there are several eigen- 
states I ] ) which have non-vanishing overlap with the 
trap node. 




l/N 



FIG. 4: (Color online) Ensemble average of the imaginary parts 
( 7 i)fl for different TV sorted in ascending order: (a) for T = 10~ 6 
and (b) for F = I. 

In Fig. [4] the average (7/} r is shown as a function of l/N 
(on double-logarithmic scales) for different TV and (a) for Y = 
1CP 6 and (b) for T = 1. One sees that the (7;)^ have a much 
richer structure than for the linear system, with evenly spaced 
nodes (compare to Fig. 1 of Ref. 11511 ). While for small TV 
a single power-law fit for intermediate values of I /TV can be 
justified (see curves for TV = 100 and TV = 200), this ceases 
to be the case for larger TV. In fact, one can identify different 



regions where different power-laws might hold, see curves for 
TV = 500 and TV = 1000. 

Due to the fact that Eq. ( TTTb is only a lower bound, relating 
exponents from power-law fits to (7; ) r to the ones obtained 
for (II(t))_R [see Fig. [3) is not very conclusive. Jensen's in- 
equality [Eq. ( TTTb l holds for the relation between (n(<))^ and 
(ji)r and thus only allows for qualitative statements. How- 
ever, it is possible to obtain quantitative relations between 
r and the average density of eigenvalues (p{^i)) b- In 
the continuum limit we get (p(j))r = £ _1 {(II(<))^}. Given 
a scaling of {H(t))R ~ one has 

(p(7))R~7' 7(JV) " 1 - (16) 
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FIG. 5: (Color online) Density of ( 7 ;)jj for TV = 100 (red circles) 
and TV = 1000 (orange triangles) and for r = 1. The solid lines 
show the scaling law which follows from the inverse Laplace trans- 
form £ -1 {n(£)} with 11(f) ~ t'^K The arrow is a guide to the 
eye pointing at the bend in p(-y) for TV = 1000. 

Figure shows the average density (/o(7;))_r (obtained by 
inversion of (ji)r) for TV = 1000 and T = 1 along with a 
power-law fit for intermediate values of 7. Indeed, the expo- 
nent of the power-law fit can be obtained from (H(t))R. As it 

turns out, the best fit is given by (p(ji))r ~ 7; 1 \ where 
ij(N) is given above for the scaling of (II(f ))# in the interme- 
diate time region. Moreover, also the slight bending of H(t) 
for TV = 1000 (see arrows in Fig.[3]l is recovered in the corre- 
sponding density /?( 7 ). 

V. CONCLUSIONS 

We have calculated the decay of the average survival proba- 
bility for coherent excitons in three dimensional topologically 
disordered systems with long-range interactions. As for a reg- 
ular one dimensional system, in the ensemble average there is 
a characteristic decay at intermediate times. We established 
a lower bound to the decay by employing Jensen's inequal- 
ity and further related the average decay rates to the average 
eigenstates of the same system without traps. Especially this 
last point might be exploited for (indirect) measurements of 
the eigenstates: For instance it is possible to control dissipa- 
tion due to the traps in Bose-Einstein condensates in optical 
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lattices [18]. Possibly also other systems with large coher- 
ence times, like optical waveguide lattices H, could be used 
to study this effect. 
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